An analytically solvable model of the effect of magnetic breakdown on 
angle-dependent magnetoresistance in a quasi-two-dimensional metal 



Andrzej NowojewskiQ Paul A. Goddard, and Stephen J. Blundelfj] 

Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom 

(Dated: February 3, 2008) 

We have developed an analytical model of angle-dependent magnetoresistance oscillations (AM- 
ROs) in a quasi-two-dimensional metal in which magnetic breakdown occurs. The model takes ac- 
count of all the contributions from quasiparticles undergoing both magnetic breakdown and Bragg 
reflection at each junction and allows extremely efficient simulation of data which can be compared 
with recent experimental results on the organic metal ft-(BEDT-TTF)2Cu(NCS)2- AMROs result- 
ing from both closed and open orbits emerge naturally at low field, and the model enables the 
transition to breakdown- AMROs with increasing field to be described in detail. 



The measurement of angle-dependent magnetoresis- 
tance oscillations (AMROs) is a powerful technique in 
the determination of details of the Fermi surfaces (FSs) 
in various reduced-dimensionality metals [TJ [21 E2 S] . In 
many cases the angle-dependence originates in correla- 
tions in the time-dependent interplanar velocity of quasi- 
particles which traverse the FS under the influence of 
the magnetic field B and hence can be efficiently sim- 
ulated by integrating up such correlations for all quasi- 
particle trajectories [3 El 12 El El DIB ID]. In high B, 
the additional effect of magnetic breakdown (MB) can 
substantially complicate this picture. This effect oc- 
curs in the FSs of quasi-two-dimensional metals such 
as that illustrated in Fig. [TTa) which is described by 
the dispersion E(k) = Tt 2 {k~ + ky)/2m* with effective 
mass m* , Fermi wave vector kp and Brillouin zone edges 
at ky == i/cp cos£. Because of the periodic potential, 
small gaps in the dispersion open up at the Brillouin 
zone edge, splitting the FS into distinct open and closed 
sections. Quasiparticles orbit around the FS with con- 
stant k z when B lies along the inter layer direction. In 
very low £?, because of Bragg reflection, only open or- 
bits [Fig. [ljb)] and small closed orbits [Fig. [ljc)] occur 
around the distinct sections of the FS. In high £?, mixing 
between the states on the two FS sections leads to MB 
at the four filled points shown in Fig. [TJa) which we term 
MB junctions. At these junctions a quasiparticle "tun- 
nels" in /c-space between the FS sections [12 , resulting 
in a single large closed orbit [Fig. [TJd)]. 

In fact for general values of the magnetic field there 
should be a superposition of all the orbits in Fig. [ljb)— 
(d) as well as many other intermediate possibilities in 
which MB occurs at some of the MB junctions and Bragg 
reflection occurs at the others. The probability p = 
exp(—Bo/B) of MB at each MB junction is parameter- 
ized by £?o, the characteristic breakdown field [12 , 13} [14] . 
For all finite, non-zero values of B (for which < p < 1) 
there is a hierarchy of complex trajectories that must be 
summed to account for all possible contributions to the 
conductivity in which MB either does or does not oc- 
cur at each MB junction. If a quasiparticle crosses TV 
MB junctions, one has to consider 2^ possible trajecto- 




FIG. 1: (a) The Fermi surface (FS) in the k x -k y plane show- 
ing the points where magnetic breakdown can occur which 
are at (k x ,k y ) = (±&f sin£, ±kp cos£) (these are called MB 
junctions). The azimuthal coordinate of a quasiparticle at 
the point labelled QP is ip. (b) An open orbit (comprising 
the repeated traversal of the /3i section), (c) Closed orbit 
(comprising the repeated traversal of a\ and a^)- (d) Break- 
down orbit (comprising a\-f3\-a2-p2)- (e) The magnetic field 
applied in a general direction leads to orbits which are on a 
cross-section perpendicular to B. Bragg reflection puts the 
quasiparticle on a different cross-section. 



ries with their correct probabilistic weightings, and this 
complicates a direct computation of AMROs since one 
has to sum over trajectories with arbitrarily long path 
lengths and hence arbitarily large values of N. In this 
Letter we describe a novel strategy to efficiently com- 
pute AMROs in a model system exhibiting MB which 
includes explicitly all these processes and we use it to 
show how different features in real data may arise. Our 
results are discussed in the context of recent experimen- 
tal work [15 on the crystalline organic metal ^-(BEDT- 
TTF) 2 Cu(NCS) 2 which demonstrated that, at high field, 
breakdown AMROs (BAMROs) could be identified in ex- 
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perimental data resulting from quasiparticles executing 
MB orbits, although until now an adequate theoretical 
description has been lacking. 

The Boltzmann equation gives the interlayer con- 
ductivity a zz = e 2 rg(E-F)(v z v z )FS as an integral over 
the FS, where v z = J °° r^e^^v^t)} dt and g(E F ) 
is the density of states at the Fermi energy. Our 
model considers the FS shown in Fig. 1(a) but in- 
cludes a very weak interlayer warping so that E(k) = 
h 2 (k 2 + k 2 )/2m* — 2t±cosk z d±, where d± is the in- 
terlayer spacing and the interlayer hopping t± is small 
(t± <C hkp/d±). For brevity, we will henceforth write 
wave vectors in units of d^ 1 and conductivity in units of 
e 2 t 2 L m*d±/h 4 7r 2 u; so that they are dimensionless. With 
B = B (sin cos 0, sin sin 0, cos 0) quasiparticle orbits lie 
in a plane perpendicular to B with angular frequency 
given by uo = cu c cos9 = eBcos9/m*. Neglecting the 
influence of t± on the quasiparticle motion (which is 
only relevant for 6 « 90°), an orbit can be described 
by k z (t) = k z + rjcos((f — (j) + | — £) where ip is the az- 
imuthal position of the quasiparticle, given at time t by 
V = + ^t, and r] = /cptan^. For later convenience, 
we measure the azimuthal angle cp anticlockwise from the 
ai-fo MB junction [as shown in Fig. [TJa)]. The interpla- 
nar velocity v z (i), which is needed to compute cj ZZi can 
be written in our units as v z (t) — sin[fe^(t)]. 

When Bragg reflection occurs for a tilted or- 
bit [see Fig. []Je)] the value of k° z jumps by 
Ak z = 2r] sin <p cos £ since only the k y value of the 
quasiparticle momentum changes and hence the quasi- 
particle continues its orbit on a different "slice" of 
the Fermi surface. We can therefore write <j zz = 
\ dk° z t d ^o sinfc(O)] /~ dipe-te-nVT sin[fc,(*)] 
where 

k z (t) = k° z + n(p)Ak z + r] cos(y> - + | - 0, (1) 

and where the term n((p)Ak z accounts for jumps in the 
value of k° z which occur during Bragg reflection and we set 
n((/?o) = 0. The k z dependence can be easily integrated 
out and we obtain 
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d(p E+((p ) / d(pE_(cp)e 



-in(ip)Ak z 



where the functions E±(x) are defined by 



E±(x) 



Q ±i?7 cos(x — </>+?r — £,)±x/ut 



(2) 



(3) 



Eq. Q yields a real expression for a zz (one can show 
straightforwardly that Im<j zz = 0). However, what 
makes Eq. ^ challenging to evaluate is that the in- 
tegrand changes depending on the path taken by the 
quasiparticle which, at each MB junction of the or- 
bit, can either undergo MB tunneling (with probability 
p = e~ B o/ B cos6>j or B ra gg reflection (with probability 
q = 1 — p): this information is encoded in the function 



n(cp) which remains constant for MB but changes by ±1 
for Bragg reflection. 

A fruitful strategy is to follow separately the motion of 
particles starting in the four different segments of the or- 
bit, only finally summing their contributions. We there- 
fore write Eq. ^ as a scalar product of vectors 



— A + • (Ai n it + x — A_), 



(4) 



where A + takes care of summing up all the initial posi- 
tions, xq handles the MB junctions, Ai n i t describes the 
initial stage of the motion up to a MB junction and A_ 
describes contributions between MB junctions. In Eq. Q 
we define 



A^ 2 

\a£ 2 

/ 



tf* dip E±(<po) 
s^ 2 ^ J^d Vo E ± (p ) 
e ^rj^ dipoE±{ipo) 



\ 



\ e ^ + 2^rJ^ dipoE±{ipo) J 



\e^»^f^ d(p E-(<p)J 



and xq is a special case of the vector 



(5) 



ft 



(6) 



representing the contribution at the MB junctions on the 
nth slice of the FS, where n is one of the integer value 
taken by the function n (<£?). 

Each component of the vectors X± and Ai n i t con- 
tributes for a specific segment of the orbit and the expo- 
nential factors multiplying some of the components are 
present in order to cancel the initial exponential damping 
of the integrand. This damping does not depend on a spe- 
cific segment but it depends on the length of trajectory 
before reaching this specific segment. This exponential 
damping is taken into account in the vector xq. This 
vector dictates the evolution of the quasiparticle's path 
and includes all the processes at the MB junctions. The 
components of the vector x n are as follows: 



n-\-l 

aqa 2 , 



(3™ = \P} e - [nAkz + bpo% + W\ 
/% = \^ e - inAkz + bpa™ + bq/% +1 , 



(T) 

(8) 
(9) 

(10) 



where a = exp(— 2^/uor) and b = exp(— (ir — 2^)/uor) are 
the increments in the damping exponential after travers- 
ing an a or /3 segment respectively. These recursive equa- 
tions encode all the information about the behaviour at 
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FIG. 2: Calculated resistivity as a function of <j) and for different values of B/Bq. In these simulations, the parameters are 
chosen to conform approximately to those appropriate for experiments on ft-(BEDT-TTF)2Cu(NCS)2(ra* = 5m e , B = 45 T, 
r = 3ps), see [151 116| . The yellow lines in (a) show the expected minima for Lebed oscillations and the blue lines show the 
expected maxima for Yamaji oscillations. 



the MB junctions. Because x n ±\ = x n e 
write them as a single vector equation 

x n = A_e- inA ^+Ex n , 



we can a 2 (p 2 — q 2 )) cos Ak z and 



(ii) 



where x n = (a™, (3™ , , 0% ) and the matrix T is given 
by a product of two matrices, one describing damping 
and the other taking into account the connection between 
orbit segments: 
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(12) 

Eq. (11) is readily solved by assigning n = to the FS 
slice initially occupied by the quasiparticle. Thus with 
X() = (I — r) _1 A_, where I is 4 x 4 identity matrix, we 
obtain 

a zz = A+ • Ai^ + A+ • £(I - E)" 1 • A_, (13) 

which is the main result of this paper. The expression 
for the r(I — r) _1 matrix is: 



/ t apr* ar*s* a 2 ps* \ 



r(i-r)- 1 



i 

N 



abps w bpr* abp 2 
2 ~ f apr 



ars a*ps t 
y bpr abp 2 abps 



(14) 



J 



where N = 1 + b 2 q 2 - a 2 (q 2 + b 2 (p 2 - q 2 ) 2 ) - 2bq{\ + 



r 
s 
t 

w 



be iAk >q 



e iAk *q 



= 1 



b(p 2 -q 2 ) 
b 2 q 2 - 2bqcosAk z - N 



= b(qe iAk 'r* + a 2 (p 2 -q 2 )s). 



(15) 
(16) 
(17) 
(18) 



The integrals in Eq. ([5| can be evaluated us- 
ing the Jacobi- Anger expansion e ^zcose _ 
Y^kL-oo i k Jk{z)z lke where J k (z) is a fe-th or- 
der Bessel function of the first kind and hence 

fE±(o)de = EZ-oo' lkJ k(±v)4 em±1/ " r) i where 

= e- ik ^-2+^)/(ik ± 1/W). This implies that 



/ J fc (±r ? )(e 2ife «e ±2 «/' 



1 \fe e ±(w-2{)/wT 
J fc ( T7? )( e 2ifc« e ±2«/^ _ 



1) \ 



1) 



V J fc (T??)((-l) fe e ± (' r - 2 «)/ WT - e 2ifc «) / 

(19) 

and similar techniques can be used to evaluate Ai n i t . Our 
results reduce to the expressions given by Yagi et al [8 
when p — > 1 or £ — > ^ while we can extract the expres- 
sions for Lebed magic angles [5] and Yamaji maxima [7] 
when p — > 0. To study the general case, we have encoded 
the solution in a computer program, separating each term 



contributing to Eq. (13) into real and imaginary parts, 



and have summed up the Bessel functions in Eq. (19), 
truncating the series at small enough J^; this involves 
typically about 200 Bessel functions in the sum. 
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FIG. 3: Calculated resistivity as a function of B/Bo for cj) — 
90° for the same parameters as in Fig. [5] 

In Fig. [2] we show the calculated resistivity as a func- 
tion of qb and for different values of B/Bo using parame- 
ters chosen to conform approximately to those appropri- 
ate for experiments on ^-(BEDT-TTF) 2 Cu(NCS) 2 . In 
these calculations B, r, and hence c^ c r, are fixed and only 
Bq is varied. When Bq is very large (small B/Bq), no 
MB occurs and the calculated p zz shown in Fig. [2|a,b) 
display a rich structure comprising features originating 
from both open and closed orbits. Some of the minima 
are due to Lebed oscillations [5] resulting from the open 
orbits [see Fig. []Jb)] and their expected angular depen- 
dence is shown in the yellow lines in the upper part of 
Fig. [2|a) and is given by tan# = mr/hp cos£ sin for 
integer n > 1 [5 . Some of the maxima are due to Yamaji 
oscillations [7J [H [17] resulting from the lens-like closed 
orbits [see Fig. [TJc)] as shown in the blue lines in the 
lower part of Fig. [2|a). The angular positions of these 
are known to be given by a "caliper" measurement of the 
FS [18] which a simple calculation can translate into the 
present geometry as 

t an6l«^— ^ x /a-sin^cosO" 1 |0-f|<£ 
J ^ k F \ (cos sin ^) 1 otherwise. 

(20) 

The values of qb at which either the Lebed or Yamaji 
oscillations dominate the calculated resistivity in Fig. [2] 
are similar to those found experimentally in ^-(BEDT- 
TTF) 2 Cu(NCS) 2 Q3]. 

As Bq decreases this structure begins to break up 
as MB starts to occur and in the limit of low Bq (see 



Fig. (2jij)) the only dominant orbits are breakdown or- 
bits (see Fig. [ljd)) leading to the observation of cj>- 
independent Yamaji oscillations with maxima at tan 6 « 
7r(n — \)/k-p. Because the breakdown probability is 
e -B /B cos e ^ ^ - g no ticeable [particularly in Figs.^h) and 
|2jj)] that p zz is more ^-independent at low than at high 
6 since MB becomes less likely as 9 increases. We note 
that at low-fields it is also possible to observe Danner- 
Chaikin-like oscillations [9] from the open sections close 
to <j> = and = § . 

The transition between the low-field and high-field be- 
havior can be studied by fixing qb and varying B/Bq and 
this is shown in Fig. [3] The series of Yamaji maxima 
(starting at « 55°), which dominates the response at 
low B/Bq due to the small lens-like orbit in Fig. [ljc) , 
give way at high B / Bo to a different series of Yamaji 
maxima (starting at 6 w 23°) due to BAMROs resulting 
from the breakdown orbit in Fig. [ljd) . The crossover 
between the two regimes begins above B/Bo w 1 which 
is when the probability of MB becomes significant. 

There is a remarkable similarity between the predic- 
tions of this model and the data of Ref. [15j [T6l It is 
expected that this approach to summing all the contribu- 
tions to the MB network model will open up new avenues 
in research on low-dimensional metals. 
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